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Abstract. We consider several tensorial wave equations, specifically the equations of Maxwell, 
Yang— Mills, and Weyl fields, posed on a curved spacetime, and we establish new energy inequal- 
C " ;) ities under certain one-sided geometric conditions. Our conditions restrict the lapse function 

, and deformation tensor of the spacetime foliation, and turn out to be a one-sided and integral 

generalization of conditions recently proposed by Klainerman and Rodnianski as providing a 
continuation criterion for Einstein's field equations of general relativity. As we observe it here 
, for the first time, one-sided conditions are sufficient to derive energy inequalities for certain 

^ I , tensorial equations, provided one takes advantage of some algebraic properties enjoyed by the 

natural energy functionals associated with the equations under consideration. Our method espe- 
cially applies to the Bel-Robinson energy for Weyl fields, and our inequalities control the growth 
of the energy in a uniform way, with implied constants depending on the one-sided geometric 
bounds, only. 
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1. Introduction 
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We consider several tensorial wave equations and, specifically, the equations satisfied by Maxwell, 
Yang-Mills and Weyl fields, posed on a curved spacetime. Given a foliation of this spacetime by 
QQ . spacelike hypersurfaces, we derive here new energy inequalities based on the observation that 

' certain one-sided bounds on the geometry of this foliation are sufficient to control the energy 

functionals associated with the tensorial equations under consideration. 

Our general purpose here is to derive estimates that depend less on norms of geometric quantities 
lO . (like the curvature) and more directly on properties such as, for example, the eigenvalue spectrum 

' of these geometric quantities. In differential geometry, one often finds that, for instance, manifolds 

with an upper curvature bound behave in a different way from manifolds with a lower curvature 
bound. Similarly, submanifolds of Riemannian manifolds with positive mean curvature are quite 
different from those with negative mean curvature. From a geometrical point of view, therefore, 
much information on the manifold is lost by simply considering the norm of a curvature quantity. 
^ , In the present paper, we identify certain minimal geometric conditions, specifically concerning 

the lapse function and deformation tensor (defined below) of the spacetime foliation; interestingly, 
these conditions turn out to be a one-sided and integral generalization of a condition proposed 
recently by Klainerman and Rodnianski as a continuation criterion for the Einstein equations of 
general relativity [9] (see also [H])- Under these conditions and by taking advantage of certain 
algebraic properties satisfied by the natural energy functionals associated with the tensorial wave 
equations under consideration, we are able to establish new energy inequalities leading to a uniform 
control of the growth of the energy, with implied constants depending on one-sided bounds of the 
lapse and deformation tensor, only. Our method especially applies to the Bel-Robinson energy 
associated with Weyl fields, which arises in the analysis of Einstein's field equations. We can 
thus identify the optimal geometric conditions that are necessary and sufficient to derive energy 
inequalities for Weyl fields on curved spacetimes. It should be recalled here that this tensor plays 
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a central role in the mathematical analysis of the Einstein equations as illustrated by works of 
Andersson [I], Andersson and Moncrief [5], Christodoulou and Klainerman [i], Klainerman and 
Nicolo [8], Klainerman and Rodnianski [9], and Reiris 

In addition, in this paper, we demonstrate that our methods may be applied to other classes of 
tensorial wave equations and, specifically, we establish new energy inequalities for Maxwell fields 
and Yang-Mills fields with compact gauge group. In contrast and somewhat surprisingly, we find 
that no such result is available for scalar fields satisfying the (minimally or conformally coupled) 
wave equation. Our results have some implications in general relativity, as we discuss at the end 
of this paper. 

Throughout this paper, we use the following notation. The manifold {Ai^^^,g) is a four- 
dimensional, globally hyperbolic spacetime, foliated by hypersurfaces TLt labelled by a time pa- 
rameter t € I, where / is an interval of the real line R. We denote by D the spacetime Levi-Civita 
connection associated with the metric g, and by R its Riemann curvature tensor. For definiteness, 
we assume that the hypersurfaces of the foliation are compact or asymptotically flat. We denote 
by N the future directed, time-like unit normal to each slice 'Ht, and by gt the induced Riemann- 
ian metric on Tit (i.e. the first fundamental form). As this is convenient, we do not distinguish, 
notationally, between the (0, 2) and (2, 0) form of the Lorentzian metric, denoting both by g. 
Using this convention, the lapse and second fundamental form of the foliation are defined by 

n:^{~g{dt,dt))-^/\ k(X,y) :=-g(DxN,r), (1.1) 

for all vector fields X, Y tangent to Tit- Next, we define the deformation tensor of the foliation by 

Tv:=CNg, (1.2) 

where C denotes the Lie derivative operator. Finally, denoting by V the Levi-Civita connection 
of the slices {Ht, gt), it is straightforward to show that 

7r(N,N)=0, 7r(N,X) = Vxlogn, 7r(X, F) = -2 k(X, F), (1.3) 

in which X, Y are arbitrary vector fields orthogonal to the normal N. Our main condition below 
will be stated in terms of tt and n, but it will be convenient to also involve k in our calculations. 

We emphasize that all solutions to the field equations under consideration are assumed to be 
sufficiently regular and have sufficient decay at spacelike infinity, so that all calculations below are 
justified. 

An outline of this paper follows. Section [2] is concerned with our main result about the Bel- 
Robinson energy of Weyl fields. In Section|31 we generalize this analysis to Maxwell and Yang-Mills 
fields and discuss the case of scalar fields. Finally, in Section 21 we discuss a class of spacetimes 
and state a conjecture for the Einstein equations. 

2. Weyl fields on one-sided bounded spacetimes 

2.1. Properties of the Bel— Robinson tensor. In all of this section (except for the final sub- 
section), we assume that the spacetime (A^'^+^,g) is Ricci-fiat, that is, in components 

Rap = 0, (2.1) 

and we analyze the Bel-Robinson energy associated to the spacetime curvature R. At the end of 
this section, we will explain how this analysis actually allows us to encompass Weyl fields on a 
given background. In (|2.ip and throughout, all Greek indices a, /3, . . . take values 0, . . . , 3 while 
Latin indices i, j, . . . take values 1,2,3. 

The following discussion depends upon the prescription of a future-oriented, time-like vector 
field which, without restriction at this stage, is taken to coincide with the normal N to the foliation 
(introduced in Section 1). In an orthonormal frame {eo, ei, e2, e^} with eo = N, whose dual basis 
is denoted by {e°, e^, e^, e'^}, the Lorentzian metric reads 

g = Vap e" ® e'^ := -e" e" -1- e^ (g) e^ -I- e^ (g) e^ -I- e^ (g) e^, 

while the volume form reads 

dVg = e" A e^ A e^ A e^ = e^^fi^s e" Ae^ Ae'^ A e*. 
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Here, as usual, we have set 

{1, ajS^S even permutation of 0123, 
— 1, af5^5 odd permutation of 0123, 
0, otherwise, 

and, for short, e^fc := eoyfc- 

The left- and right-hand Hodge duals of the curvature tensor R, by definition, are 

respectively. We recall that if the metric g is Einstein (i.e. the Ricci tensor Ric is proportional 
to the metric g), then "^R and R* coincide. In addition, when g is Ricci-fiat as we assume in this 
section (see (|2.ip ). then the tensors *R and R* possess the full symmetries of the curvature tensor 
R. (For details, see for instance JD, Sec. 4.6].) 

From the curvature tensor and its dual, we may then construct its electric and magnetic parts 
with respect to the vector field N, i.e. 

E(X,y) := g(R(X,N)N,r), H(X,y) := g(*R(X, N)N, F) , 

respectively. By the (skew) symmetries of the curvature tensor R it follows that E and H are 
symmetric and tangent to the foliation hypersurfaces T-Lt, i.e. 

E(N,.) = 0, H(N,.) = 0. 

Moreover, since the metric is Ricci-fiat, both E and H are trace- free, in the sense that (5'-'E(ei, e^) — 
5*J"H(e„ej) =0. 

To derive certain relations between the quantities R, *R, E, and H, it is convenient to define 
the reference Riemannian metric, 

gN :=g + 2g(N,.)®g(N,.), (2.2) 

associated with the Lorentzian metric and the given vector field. This Riemannian metric deter- 
mines the corresponding norms | • and inner products ( • , • )gj^ on tensor bundles over Ai^^^. 

Proposition 2.1 (Decomposition of the Riemann tensor). Relative to an orthonormal frame 
{eg, ei, 62, ea} with Gq = N, one has 

^ijkl — eijm^kln-^mm ^ijkl — *^ijm^kln-^mn- 

Moreover, with respect to the Lorentzian metric g and the Riemannian metric gN one obtains 

|2 

■IgN 

respectively. 



Rp = 8 (|Ep - |Hp), |R|2^ ^ 8 (|Ep + |Hp), (2.4) 



A proof of this result is provided shortly below. To the Riemann curvature tensor, we can 
associate a notion of energy, i.e. the Bel-Robinson tensor 

Qal3fS ■= Q[R]q^7i5 := RaXjfj.Rl3^S^ + * RaXffi* Rfi^ ■ 

It is straightforward to check that this tensor field is totally symmetric and trace-free. Moreover, 
the components of Q[R] are related to the electric and magnetic parts of the curvature, as observed 
in 131 15 . These relations are conveniently expressed by setting 

A-B^A^pB''^, {AAB),,=eJ'<Ap^Bsj, 

and 

{A X B)afi =e^'^^ep^"'A^^Bsu + B)g^p - i(tr A)(trB)g„^ 

^A^^B^p + A^T^^c B)gc,0 + ^(tr A)(trB)g„/3 - (tr A)B„^ - {trB)A^p. 

Observe that the x operation is symmetric. 
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Proposition 2.2 (Decomposition of the Bel-Robinson tensor). Relative to an orthonormal frame 
{60,61,62,63} with 60 = N, the relevant components of the Bel-Robinson tensor take the form 

Qoooo = |Ep + |Hp, Qiooo — 2(E A H)^, 

1 (2-5) 
g,,oo = -(E X E),, - (H X H),, + -(|E|2 + \U\^) g,^. 

Proof of Provositions [2J\ and \2.Sl To derive (|2.3p we observe that, by definition, 

* RijkO = -^f-ij^^ Rpiuka 

— ^^j RomkO — ( ^ijm^{ RmOko) — ^ijm-Ejjik- 

The property = — id implies 

RijkO = —* i* Rijko) = {^ijmRmOko) 
^ijm RmOkO ^ijm-^mk- 

The identities in the first line of (|2.3p then imply those in the second line, more precisely 

^ijkl — 2 ^^fkl ^ijm^Omkl 

— ^ijmRklraO — ^ijm^klnHmn 

and 

Rijkl = ^* {* Rijkl) = ^* { — ^ijmRomkl) 
^ijm RklmO ^ijm^kln-^mn • 

Due to the symmetries of the Riemann tensor, the above quantities suffice to derive 

Rp — Rap^sR"^''^ = '^RiojoR^^'^^ + ^RijkoR^'^'^^ + RijkiR^'''^^ 

— ^-Ej^j E ^ ^^ijra^ijnHkmH ' ^ ^ijm^kln^ijo^klpEmnE ^, 

thus 

= 8f|E|2- IHPI 



On the other hand, to handle IR-lg,^, we observe that, in the calculation above, we do not get a 
factor of —1 if we lower the index on the second term, and this then changes the sign of the 
term |Hp. This completes the derivation of (|2.4p . 

Next, considering the Bel-Robinson tensor Qap-ySj we immediately deduce from the definition 
that 

Qoooo = RoiOjRo^o'' + RoiOj Ro^ 0'' 

= E,jE'^ + HijW^ = |E|2 + |H|2, 
while all the other terms can be computed in a similar fashion. □ 

2.2. Bel Robinson energy inequality under a one-sided bound. On a Ricci-fiat spacetime, 
the Bel-Robinson tensor has the divergence-free property 

Ti'^QcMS = 0, (2.6) 

as observed in Penrose and Rindler [TOl Sec. 4.10] and Christodoulou and Klainerman [4, Propo- 
sition 7.1.1]. This property suggests that by integration one should be able to control the Bel- 
Robinson energy on an arbitrary slice of the foliation in terms of its values on some "initial" slice. 
More precisely, we introduce the total Bel-Robinson energy at time t as 



Q[R]n, ■■= [ Q[R](N,N,N,N)dV^g 

JHt 



and we now derive a uniform estimate for Q-Ut which solely involves a one-sided bound on the 
deformation tensor of the foliation. 
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First, we fix some notation and, by adopting an orthonormal frame with Gq = N, we express 
the deformation tensor defined in (|l.ip - (ll.3p in the form 



Vn 

Vn -2nk 



(2.7) 



We also define the auxihary tensor field 

A(nk) (tr(nk)) ids -2nk, 

and the cubic polynomial of the real variable A 

P„^(A) := det (A idg -A(?ik)) + (A(nk) - A ids ) (Vn, Vn) , (2.8) 

referred to as the critical polynomial of the spacetime foliation. This polynomial will arise naturally 
from the expression of the Bel-Robinson tensor; see (|2.16p . below. 

Definition 2.3. The largest real root of the polynomial Pm^ is called the critical root of mr and 
denoted by p{mv). 

Theorem 2.4 (Bel-Robinson energy inequality under a one-sided bound). Given any vacuum 
Einstein spacetime endowed with a foliation {'Ht)t£i with lapse function n and deformation tensor 
IT, one has 

for all ti,t2 G /, where 

t2 



Kn-7v{ti,t2) ■■= / supp{mr)dt 
Ju Ht 



Ui -Ht 

For the proof of this result, we will derive (in Section [2. 3[ below) the algebraic inequality 

- ^ Qapoo mr"^ < p{nn) Qoooo- (2.9) 

At this stage, we only check that (|2.9p implies the energy inequality stated in Theorem 12.41 
Namely, fixing any ti,t2 G I such that ti < t2 and applying Stokes theorem to the vector field 
Pq — Q\R]aP'ysNl^ on the manifold with boundary A^fti.ta] := UtG[ti t2] obtain 



dV^ 



which implies 



1/ Qc,0^s7T''^'N^N'dV^= f D"P„ 
<QHt^+3 / p{nn) Qqoqo dVg^ds 

Jti JHs 

/■*2 

< Qut, + 3 / sup(p(n7r)) Qn^ds, 

where, in the first inequality, we have used (|2.9p . Since sufficient regularity has been assumed on 
all solutions under consideration, we deduce that the derivative of Q-Ht a.t t = ti satisfies 



lini < 3sup(p(n7r)) Q-H 

12— >tl to — t-\ 



t=ti '^^'1 ''2 ~ '•1 Wt 

Since this inequality holds for arbitrary ti € I, we may integrate it and arrive at the energy 
estimate of interest: 

Q-Ht^ < Q-Ht, exp supp(n7r) dt^ , ti < tz- 
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2.3. Derivation of the optimal condition. This section is devoted to deriving the algebraic 
inequality (|2.9p and, more precisely, 

- I Qa^m mr"^ <CQoooQ (2.10) 

for some constant C > depending on certain assumptions stated below. By using the algebraic 
properties on Qai3oo stated in Proposition l2.2[ above, we now derive conditions (on the eigenvalues 
of mr) implying ()2.10|) . We discuss the general case first, and then consider special cases of 
independent interest. 

Proposition 2.5 (Main algebraic condition on tt). The product — ^Qq^oo ""■"^j decom- 
posed as in (|2.7p . can be controlled by the (double) eigenvalues ai < 02 < 03 of the 6 x 6-matrix 
n, given by 

/A(k) ~S{7t)\ ^ 



n^="(,5(7r) A(k) j' S,,i7v) -.^ e,,,7: 



where the 3 x 3-matrices A and S are symmetric and skew- symmetric, respectively. More precisely, 
the estimate (j2.10p holds with C = 03, and the eigenvalues ai < a2 < 03 of II are nothing but the 
roots of the polynomial Pn-w (defined in ()2.8p ). so that 

03 = p{mT). 

Proof. From Proposition 12.21 we compute 



2 



so 



(% (|E|2 + |H|2) - 2 {E,kE,k + H,kH,k)) k'' - 2e,jkEjiHikTr"' 
trk (|E|2 + |H|2) - 2tr (k (E^ + tf)) - 2S,kE,iHiu, 



-^Oo^oot"'^ = tr (A (e2 + H2) - 25HE) 

= tr (EAE + HAH - E^H + HS-E) (2.11) 



and, equivalently. 



'^Qa0oon7r"^=ntr('(E,H)('^ /V^ 



tr( (E,H)n(^ 



(2.12) 



which provides us with a rather explicit expression for the left-hand side of (|2.10p . 

Observe that the 3x3 matrices A and S are symmetric and skew-symmetric, respectively. The 
linear map 11 given in the statement of the proposition is a symmetric 6x6 matrix. Moreover, 



if we have an eigenvector of 11 with eigenvalue a and write it in the form where u,v £ 



p3 



V 



-u 



then it is obvious that is a distinct eigenvector of 11 associated with the same eigenvalue a. 



Therefore, the matrix 11 has (at most) three distinct and real eigenvalues, each of which appears 
with multiplicity two. Denoting these eigenvalues by ai < 02 < 03, from (|2.12[) and p.Sp . we 
deduce: 

- ^(3Q/3oo^^7I■"^ < max(ai, 02,03) tr ^(E H) ^ = 03 Qqooo- (2-13) 

Next, we need to study the eigenvalues of the matrix 11 in terms of tt. First, we consider the 
case that the matrix A is diagonal, so of the form diag[^, B, C], and we set 

a := -5*23 = Vi logn, 6 := -53i = V2 log n, c := -512 = V3 log n. 
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We therefore have 



n 














c 


-b\ 





B 





— c 





a 








C 


b 










— c 


b 


A 








c 





—a 





B 







a 











c) 



(2.14) 



and a simple calculation leads us to 

det(n- A ids) = P(A)^ (2.15) 

where 

P(A) :=A3 - X^n{A + B + C) + \n^{AB + BC + CA- {a^ + b^ + c^)) 



{^ABC + Aa'^ + 56^ + Cc^) 



(2.16) 



= det (A ids -nA) + (n^ - A) ii^a^ + (nS - A) n^b^ + (nC - A) n^c^ 

= det (A ids --"A) + (nA - A ids) ( Vn, Vn) . 

Hence, the eigenvalues of the matrix 11 are given by the roots of the cubic polynomial, Pn-K 
defined in (|2.8p . From the fact that the latter expression in (|2.16p is S0(3)-invariant, we deduce 
that the same conclusion holds even when A is not diagonal. The discussion above implies that 
the polynomial Pn-^ has three real roots ai < 02 < 03, and the result then follows from (j2.13p . □ 

Remark 2.6. In ()2.16p . only the squares of the coefficients a, b, c appear. Therefore, to bound the 
roots of the polynomial Pmr, we will need to impose a bound on |Vn|. However, the fact that the 
polynomial Putt (A) is invariant under the change of sign of a, 6, c will later work at our advantage 
and allow us to carry over our analysis to the case of Maxwell and Yang-Mills fields with only 
minor changes. 

2.4. Two special cases. Although it is possible to write down the roots of a cubic polynomial 
explicitly, in most cases it is more useful to estimate these roots by simpler quantities. We first 
consider two illustrative special cases of vanishing k and vanishing Vn, respectively. 

Corollary 2.7 (The case of vanishing k). Assume that the second fundamental form of the foli- 
ation slices vanish identically and, therefore, tv is of the form 

V log n 
V log n 

Then, the condition (|2.10p is satisfied with C = |Vn|. 



Under the assumptions in Corollarv l2.7l we obtain a stronger conclusion and, actually, both an 
upper and a lower bound for the Bel-Robinson energy. (Cf. the two inequalities (|2.17p . below.) 

Proof. By setting A and Sjk £ijk^i log?i in Proposition 12.51 '^^ obtain 



Qa/^oo"-7r"'^ = "^tr ^(E,H) ) ) < max(ai, 02, 03) Q, 

The eigenvalues ai , 02 , 03 of the symmetric matrix 

'0 -S' 



0000- 



n 



S 



are the roots of the characteristic polynomial (I2.16p . i.e. ai — — |Vn|, a2 = 0, and a^ — |Vn|. 
Therefore, we find 

- I V?i| Qoooo < - ^Qa/soo ^^7^"^ < | Vn| Qoooo, (2-17) 

and Qajsoonn"^ is controlled by |V7i| and Qoooo- D 
From now on, we denote by 

fci < fca < fca (2.18) 
the eigenvalues of the second fundamental form k of the foliation slices. 
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Corollary 2.8 (The case of spatially constant lapse). Assume that the lapse of the foliation is 
constant and, therefore, tt is of the form 




Then, if the inequality 

n{k2 + h) <C + nki (2.19) 
is satisfied for some real C , then the condition (|2.10p holds with this constant. 

The condition (|2.19p implies that 

nki < nk2 < nk^ < C, 

so that the principal curvatures of the hypersurface (weighted by the lapse n) are bounded 
above. Geometrically speaking, this condition means that the normal geodesies to the hypersurface 
T-Lt are not focusing inwards too quickly. It seems reasonable that such a condition is required if 
we wish for our foliation not to collapse. 



Proof. Since the second fundamental form is (real-valued and) symmetric, we may, without loss 
of generality, assume that k is diagonal of the form diag[fci, fc2, fca] with, say, ki < fc2 < ^3. We 
then deduce that A(k) = diag[Ai, A2, A3] with 

Ai -fci + ^2 + fcs, A2 := fci - fc2 + fca, A3 fci + fc2 - ^3- 

Observe that A3 < A2 < Ai and, in view of equation (|2.8|) , that the eigenvalues of the matrix 11 
are n\i. As such, the condition (|2.10p is satisfied if and only if 

nAi < C, 

which gives (|2.19p . □ 



Remark 2.9. Another way to view our result is to relax the ordering condition on the eigenvalues 
and impose that Ki :~ nki — C < 0. Defining the vectors 




then the condition (|2.19|) may be rewritten in the form (K, ii)^^ < 0. Therefore, the vector K 
lies in the intersection of three half-planes through the origin in M'^. As such, nk lies in a non- 
compact, triangular cone in with vertex at the point (C, C, C). Although the eigenvalues of 
nk are bounded above, they may all be arbitrarily large as long as they remain within this cone. 
Observe that this cone contains the case where nk is a multiple (< C) of the identity matrix 
and, therefore, contains matrices that are not too far from a multiple of the identity matrix. The 
bound under consideration is very different from the bounds considered in [9l [12] which impose a 
restriction on |k| and, hence, impose both upper and lower bounds on the principal curvatures of 
the hypersurfaces of the foliation. 

Remark 2.10. Under the conditions in Corollarv 12.81 assume also that the foliation consists of 
slices with constant mean curvature, that is, trk = < on for all t G I. If k satisfies (|2.19p . then 
2nki{t) + C > ntrk^ nt and, therefore, 

nki{t) >^{nt~C). 

It follows that, as long as the lapse is bounded away from zero, then k is automatically also 
bounded from below as well as from above. 
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2.5. Derivation of sufficient conditions. For general tensors tt, we may combine the results 
of Corollary 12.71 and 12.81 and arrive at a simplified condition on tt, that is, the condition (|2.20p . 
below, which is in accordance with (|2.21[) derived in the general case. 

Proposition 2.11 (A sufficient condition for general foliations. I). A sufficient condition for ()2.10|) 
to hold is that 

n{-ki + k2 + k3) + \Vn\<C, (2.20) 
which is automatically satisfied under the stronger condition 

C 

max (n (— fci + fc2 + k^) , | Vn|) < — . 

Proof. From (|2.1ip . we have 

Qa/3oo7r"'^ = n tr (EAE + HAH) + n tr [HS'E - E5H) . 

We now estimate the two terms on the right-hand-side of this equation separately. CoroUarv 12.81 
implies that the first term is bounded above by n (— fci + k2 + k^) Qooooi while Corollarv l2 . Tl implies 
that the second term is bounded above by | VnlQoooo- D 

Using estimates for eigenvalues of symmetric matrices, we deduce the following assumptions 
that are sufficient to bound Qapoo'^"^ without having to compute all eigenvalues explicitly. 



Proposition 2.12 (A sufficient condition for general foliations. II). By defining 

r2"|+n|7r30|, 

.,301 I ^|„10| 



ai n (-fci + k2 + k^) + n|7r^"| + n\Tr^ 



n (fci - fc2 + fcs) + n|7r-^"| + n|7r^"|, 
n(fci + fc2-fc3) + n|7ri0|+n|7r20|, 



0'2 
f73 

the condition (j2.10p holds with 

C = max((Ti, (72, (T3). (2-21) 
In particular, this holds if nk. satisfies (|2.19p and |Vn| is bounded. 

Proof. Gershgorin's circle theorem allows us to bound the spectrum of the symmetric matrix 11 
and, precisely, states that the eigenvalues of 11 (with the notation (|2.14p ) are contained in the 
(closed) balls B„|f,|+„|c|(nA), Bn\a\+n\c\{n-B) and Bn\a\+n\b\{nC). Thus, an upper bound on the 
eigenvalues is given by max(cri, cr2, cts)- Together with Proposition 12.51 this leads to the desired 
estimate ([TTO)) . □ 

2.6. Weyl fields on curved spacetimes. Our previous results are now generalized to Weyl fields 
defined on a fixed Lorentzian manifold (A4"^+^,g). Recall from Christodoulou and Klainerman [H 
Chapter 7] that a Weyl Geld is a (0,4)-tensor field, say Wap^s, that has the same symmetries as 
the Ricmann tensor and, in addition, is trace-free, i.e. 

Wpo.S = 0. 

Again, the left and right duals, 

are equal, and are also Weyl fields. Finally, we define the Bel-Robinson tensor of the Weyl field 
to be 

Q := Q[WU^s Wc^x^f^Wp^^ + *Wa,xjt.*W p^s'' ■ 
We assume that the spacetime A4^^-^ is endowed with a foliation as stated in the introduction, 
and we define the totai Bel-Robinson energy of the Weyl field at time t to be 



Q[W]n,:= [ Q[W](N,N,N,N)dV^g,. 

JHt 



10 



A.Y. BURTSCHER, J.D.E. GRANT, AND P.G. LeFLOCH 



In the special case that the Bel-Robinson tensor associated with W is divergence-free, we imme 
diately obtain an energy inequality for Q[W] along the same lines of the proof of Theorem 
specifically we obtain 



2- 



Observe that we no longer require that the background metric is Ricci-flat. However, in general, 
div(Q[W]) 7^ and, instead, we have the following conclusion. 

Theorem 2.13 (Energy inequality for Weyl fields under a one-sided bound). Let {Ai^^^, g) be a 
spacetime endowed with a foliation {'Ht)tei with lapse function n and deformation tensor tt and, 
for all ti,t2 G I with ti < t2, define 

Kn7r{ti,t2) -.^ / supp{mv)dt, 

where p(rnv) is the critical root of Pn-K (cf Definition \2.3\) . Then, any Weyl field W defined on 
A4^~^^ satisfies the energy inequality 

Q[W]„,^ <Q[W]„,^e3^'"'(*-*^) - / ' e3^-(*^*-) / n {div Q [W])ooo dV^.dt, t^ < t^. 

Jti JHt 

Proof. Recalling that div((3[W]) ^ in general, we find the energy identity 
Q[W]«,^ = Q[W]n,^ -I r I nQ^p^o^^f" dV^Jt 

-T I n(divQ[W])ooodl^g,di (2.22) 

for all ti,<2 G I with ti < <2- Arguing as in the passage after Theorem 12 .41 we deduce from the 
inequality (j2.22p that the derivative of Q at ii e / satisfies 

and, therefore, for all ti < t 



<3sup(p(n7r))Q[W]„,^ - / n (div Q[W])ooo d^g,, 
J (e-3^""(*-*)Q[W]„,) < -e'^K^M j „ (div Q[W])ooo d^g,- 



Integrating in ^ G [ti , ^2] , we deduce that 

e-3^""(*^'*^)Q[W]„,^ -Q[W]„,^ <- / ' e-3^""(*i^*) / n (div Q[W])ooo rfVg.di, 

as required. □ 

Remark 2.14. We cannot proceed any further without placing additional restrictions on the 
divergence of the Bel-Robinson tensor of the Weyl field, a quantity which, in the context of the 
Einstein equations, is determined by the matter content of the spacetime under consideration. Such 
terms with non- vanishing divergence do arise, even in the context of the vacuum Einstein equations, 
as error terms generated by commuting Bianchi equations with vector fields (as developed in [4]). 

3. Maxwell and Yang-Mills fields on one-sided bounded spacetimes 

3.1. Maxwell equations. We now consider a vacuum Maxwell field, represented by its curvature 
(or field strength) F, defined on a four-dimensional Lorentzian manifold {M.^'^^,g,) foliated as 
before. By definition, F satisfies the vacuum Maxwell equations 

dF = 0, d*F = 0. (3.1) 

The energy-momentum tensor of the Maxwell field takes the form 

Tap — Fa^Ffj'^ — -g^^lFlg. 
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Letting e" be an orthonormal coframe, with dual frame Bq,, where we take Gq — N, then we divide 
the components of F into electric and magnetic parts: 

We define the total energy of the Maxwell Geld at time t to be 

Q^W[F]„, := / T(N,N)dyg, = 1 / {\Ef + IB^ dV^,. 

A straightforward calculation in the above orthonormal frame then shows that 

(T, tt) - i trTT (|E|2 + |B|2) - ^^J- {E,E, + B,B,) - 2^°^ (E x B), 

Letting Sij := eijk'^^'^ and noting that Tr^ = ~2kij, we deduce that 

(T, tt) = - trk (|E|2 + |B|2) + 2k'J {E,Ej + B,Bj) - 2SjkEjBk 
= tr [(E^ + B^) (2k - tr k idg)] - 2 tr [ESB] . 
As before, letting A(k) := (trk) ids —2k, and noting that S is skew-symmetric, we deduce that 



-i(T, tt) = i tr (EAE + BAB + ESB - Bfi-E) . 



We therefore have 



(i(E,B)(^, f)(^ 



which leads to similar formulas to the ones found with the Bel-Robinson tensor. 
We therefore define 

' A 5^ 



n.-,. A 

Note, from the explicit form of the characteristic polynomial of 11 given in (|2.15p and (|2.16p . that 
changing the sign of the matrix S does not affect the eigenvalues of the matrix 11. It follows that 
the eigenvalues of the matrix 11 coincide with those of 11. We therefore deduce that 

-in(T,7r) < imax(a,)tr (^(E,B) =max(a,)Too. 

Following the same proof as in our estimate for the Bel-Robinson energy, we therefore deduce the 
following result. 

Theorem 3.1 (Energy inequality for Maxwell fields under a one-sided bound). Let (A^'^+^, g) he 
a spacetime endowed with a foliation {'Ht)tei with lapse Junction n and deformation tensor tv and, 
for all ti,t2 I with ti < t2, define 

Kmv{ti,t2) ■■= / supp(n7r) dt, 

where p{mT) is the critical root of Pn-K (cf. Definition \2.S\) . Then, any vacuum Maxwell field 
defined on JH^^^ satisfies the energy inequality 

3.2. Yang— Mills equations. As should be clear, the calculations for the Maxwell field may easily 
be generalized to the case of Yang-Mills fields F with compact gauge group, G. In this case, the 
field strength is a two- form taking values in the associated Lie algebra g. Since G is compact, 
the Cartan-Killing form k is negative definite. Letting T* (i = 1, . . . , dimG) be an orthonormal 
basis for g with respect to the Cartan-Killing form, then the energy-momentum tensor of the 
Yang-Mills field takes the form 



dimG 
i=l 



Faj'Fp''' - lg„^F^5^F^*A . (3.2) 
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The electric and magnetic parts of the curvature, along with the Yang-Mills energy are defined 
in precise parallel with the case of Maxwell fields. As for Maxwell fields, we find that 

4»(T..)=i./|f.,((E-,B-)(_'\, 

From this result, again following the proof from the Bel-Robinson energy estimate and denoting 
by Q^'^[F]-Ht the total energy of the Yang-Mills field F at time t, we deduce the following result. 

Theorem 3.2 (Energy inequality for Yang-Mills fields under a one-sided bound). Let (A^'^+^,g) 
be a spacetime endowed with a foliation {'Ht)tei with lapse function n and deformation tensor tt 
and, for all ti,t2 G I with ti < t2, define 

rt2 



Kmv{ti,t2) ■■= / supp{mT)dt, 



where p{mv) is the critical root of Pmr (cf. Definition \2.3\) . Then, any Yang-Mills field (with 
compact gauge group) F defined on A4^^^ satisfies the energy inequality 

Q^'^iFjn., < Q^'^mn,, e^""(*-*^), t, < t^. 

3.3. Scalar wave equation. We now discuss the scalar wave equation, with the intent of also 
sharpening the standard energy inequality. As we shall discover, however, this equation does not 
have a "sufficiently rich" algebraic structure for an analogue of Theorem 12.41 to hold. 

As before, (A^'^^^ , g) is a Lorentzian manifold, but we allow it to be of arbitrary dimension d-l- 1 , 
and we assume a foliation by level sets T-it of a time-function, as described in the introduction. 
Let ip: M^^^ ^ R be a real- valued scalar field on A^'^+^ that satisfies the wave equation 

□g(^ = 0. (3.3) 

Recall here that all solutions under consideration are smooth and have sufficient decay at spatial 
infinity. The energy-momentum tensor of the scalar field Lp reads 

T := (g) diyj - i {dip, dip) g, 
and the total energy of the scalar field at time t is 

Q^'^Mw. T(N,N)dFg, = i / \T>p\l^dV^,. 

JUt ^ JUt 

We then have for ti < t 

S^^M«. - Q^^M«., = - / Co (T-'^iV^) dV^ = -\f (T, 7r)g dV^, 

where (T, 7r)g denotes the inner product (of the (0, 2) tensor fields) with respect to the Lorentzian 
metric g. 

Let us first recall the standard inequality. We would normally estimate the second integral by 



/ (T,7r)dl/g< / |T|g,|7r|g,dFg 



and, since |T|gj^ < |D(^|g deduce that 



/ (T,7r)dFg < / \T>p>{t)\lJ7,{t)\^^dV^^ I I \T>p\lJ7,\^^ndV^ds 
TsupdnTTl) / \Dp{s)\l^ds^2 [\xx^{\m,\)Q'^''[ip]n^ds. 



< 

Gronwall's inequality leads us to the standard conclusion 



Q'^^Mw. <exp (2 rsup|n7r|gds)QWE[^]^^^ 
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Turning now to the derivation of a possibly weaker restriction on the deformation tensor, we 
need to consider the quantity (T, 7r)g more carefully. We have 

(T,7r)g = T^pir'^^ = tt"'' (^D^^D^^ - ^5a/3|^¥'lg) 

= ~\ (t'^g ^) (gradg V, gradg ^p), 
where we have defined the vector field gradg by, for arbitrary vector fields X, 

(X, gradg <^)g =X((p). 

Let L be the operator defined by 

(X, Ly)g, =n{^ (trg tt) g - tt) (X, F), 

or, equivalently, in local coordinates, 

L^fi ■= ng^ Q (trg tt) g^^ - -k^^ . 

In terms of an orthonormal frame ea with eo := N, L has components 
I/°o = --^(trgTr), L*o = -nTTio, 



-niTio, Uj = n ^^(trg TT)Sij - nij^ 



and, in particular, tr L = ^{d — 3) trg tt. 

We denote the eigenvalues of L by Aq, . . . , Ad, and then deduce that 



1 , 1 

2 
1 



< -max(A„)|gradg93|gj^ =max(Aa)T( 



00- 



^ ^. .^j,^ I ^ Lg -^ig 

Our energy inequality for the wave equation = therefore becomes 

< (^J max(Aa)Toodyg,^ ds 

< / max(max(AQ(s))) I / Tqq^s) dVg\ ds 
Jti \JHs J 

= (^max(max(Aa(s)))^ Q^^M^.rfs. 

As such, if we impose that the eigenvalues of the operator L are uniformly bounded above, in 
the sense that 

{X,LX)^^<C\X\l^, 

then we have the inequality 
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Since Q^^l'pl-Us 0, it follows that we have the energy inequality 

Q^'^Mh, < Q^^M«., exp (^l* C{r)dr 

As the following example shows, however, the condition that the operator L is bounded above 
implies that the second fundamental form is bounded both above and below. Consider the case 
where the lapse, n, is constant on the hypcrsurfaces Ht, i.e. n = n(t). The only non- vanishing 
components of the deformation tensor are nij — — 2fcy . We then have 

-in(T,7r) = n ( k'^ -UrksA D^^Dj^ + ^ trk{Do^f 



Let A'^ := k^^ - ^trkS'^. Observe that A is a symmetric, real {d x d) matrix, which may, 
without loss of generality, be assumed diagonal of the form A = diag[Ai, . . . ,Xd]- Observe also 



that tr k = ^7jr2 ^^'^ then have 



-^{T,tt) ^ntr{AD^(E)Dip) - -^tr A {Do^pf 

d 

= J2K{D,^f-^trA{Doipf 

i=l 

= {Do(p, Ddip) A (Dqlp, Dd^f , 



where 



A :— diae 



d~2 '^i'---'^" 



In order to bound — in(T, tt) above in terms of Tqo, we need the eigenvalues of the matrix A to 
be bounded above. This implies that we require the existence of a constant C such that 

Ai,...,Ad<C, (3.4) 

but also 

Ai + •• ■ + Ad > -C(d-2). (3.5) 

However, the bounds (|3.4[) and (|3.5p together imply that 

Ai,...,Ad> -C(2d-3). 

It follows that, since A is bounded both above and below, that the second fundamental form of 
the foliation, k, must also be bounded both above and below. 

As such, our method does not generalize to the scalar field when we only have a one-sided bound 
on the deformation tensor. In particular, our methods depend strongly on the algebraic structure 
of the energy momentum tensor. The trace-free nature of the Bel-Robinson tensor and the energy- 
momentum tensor for Maxwell and Yang-Mills fields suggests that one might fare better with a 
conformally coupled scalar field, which satisfies 

□ g<y9 + C(d)i?gV3 = 0, 

where C{d) = and i?g denotes the scalar curvature. The energy momentum tensor of such a 
field takes the form 

T^p = (1 + 2C{d)) Da^Db^ - ig„^ (1 + AC{d)) \T>ip\l 

+ ^gapC(d)Rgip^ - C{d)Raf!f^ + 2C{d)ipD^Dp'p - 2C{d)ga0'fU^'p. 

Unfortunately, it is now not possible to control — n(T,7r) due to the terms of the form ipDaDpip 
in the energy-momentum tensor. 
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4. Example and conjecture 

4.1. Example of spacetimes. Our estimate in Theorem 12 .41 shows that, mider certain one-sided 
bounds, the Bel-Robinson energy of a vacuum solution of the Einstein equations cannot grow too 
quickly. If one were to assume a two-sided bound on the deformation tensor, then one would find, 
by a similar argument, a lower bound on the Bel-Robinson energy, showing that the energy, in 
addition, cannot decrease too quickly under a two-sided bound. In this section, we show that one 
can always choose a Lorentzian metric that satisfies a one-sided bound of our type, but for which 
the Bel-Robinson energy dies off as quickly as desired. More precisely, we have the following 
result. 

Theorem 4.1. Given arbitrary e > and A G (0, 1), there exists a foliated, four- dimensional, 
Lorentzian manifold {A4^^^,g), with A4^^^ = IJtG/^*' with the property that for to G I 

Hence, the Bel~Robinson energy dies off by an arbitrarily small factor (i.e. X) in an arbitrarily 
short time (i.e. e). 

Proof. Our manifold is based on the Kasner metric. Recall that the Kasner metric takes the form 

g = -dt^ + t^'^dx^ + t^^dy^ + t^'^dz^, 

where we take t > Q and (x, y, z) to be coordinates on a three-dimensional flat torus, T^, with the 
property that ^j,^ dxdy dz = 1. We choose to take a < /3 < 7. The above metric is then Ricci flat 
if the conditions 

are satisfied. Excluding the case where a — P — Q, ^ — 1 which gives a flat metric, then we 
necessarily have a < and 7 > i. A straightforward calculation shows that the eigenvalues of 
the second fundamental form of the surface St (of constant t) are —a/t, —ji/t, —^/t. As such, the 
first root is unbounded above as i — ^ 0, while the last is unbounded below. It turns out that the 
behaviour of the root — 7/t will be the most relevant to our discussion. 

Noting that the lapse of this metric is equal to 1, the operator A(k) defined earlier has eigen- 
values 

Ai :- A2 :- A3 :- 

It follows that the critical root is 

p(mT) = — ^ — . 

A calculation of the curvature then shows that the magnetic part of the Weyl tensor of the Kasner 
metric vanishes, while the electric part is given by 

t^ |Ep ^a^{a^ If + - If + 7^(7 - if. 

It follows that the Bel-Robinson energy is 

Qn, = a\a - if + - if + 7^(7 - l)^ t > 0. 

Our metrics are constructed as follows. Let e, A > be given, with A < 1, and define the 
quantity 

eAi/3 

"^="T^AV^<"- 
In terms of a and for all t > a, we define the translated metric 

ga -dt^ + it~ af" dx^ + {t~ af^ dy^ + (t ^ af" dz\ 
Then, the Bel-Robinson energy of the set t = constant > a is given by 

_ aHa-lf+(3H/3-lf+^^i^~lf 
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In particular, letting tQ = 0, we find 



^Vl"^— V' = A 

and it follows that the metric ga, with = 0, has the required property. □ 

Remark 4.2. The metric has constant lapse, so the critical root of the deformation tensor is 
given by 

( A 27-1 

p{mr) = . 

t — a 

For t > to > 0, we therefore deduce that 

* t 
p{n7z){s) ds — (27 — 1) log 

to 



to - a 

This is bounded for any finite t > to, and therefore our inequality in Theorem \2.4\ applies on any 
finite t-interval. 

4.2. Einstein equations. This paper stems from recent work on the vacuum Einstein equations 
by Klainerman, Rodnianski, and Wang [9j 112) . in which breakdown criteria were proposed for 
regular solutions to the Einstein equations. Building upon the strategy already adopted by Grant 
and LeFloch [Sj |6] to deal with the injectivity radius of Lorentzian manifolds, we required in the 
present paper solely one-sided geometric bounds on the spacetime. Our method should in principle 
be applicable to the Einstein equations. 

Recall that in harmonic coordinates x^ , satisfying by definition g'^^T'^^ — where T'^^ denote 
the Christoffel symbols, the vacuum Einstein equations take the form of a tensorial system of 
nonlinear wave equations 



^ g'^^d^dsgap + Fapig, dg). 



= i?a/3 = 2 

The metric (/q^ = ma/3 + Ua/i is typically expressed as a small perturbation u of the Minkowski 
metric m, and the reduced vacuum Einstein equations above are rewritten as a nonlinear wave 
equation for u, i.e. 

□,„u = A^(u, du,d^u), 

where Dm := rn^^dadjs, and the nonlinearity N consists of terms of the form F(u) • u • d^u and 
F(u) • du ■ du. Therefore, relative to wave coordinates, breakdown of smooth solutions can not 
occur on some interval [to,ti) (see [7J Sec. 6.3] and [5J Sec. 3]) as long as 

[ \\dg{t)\\L^(^-Ut)dt < 00. 

Jto 

On the other hand, by taking advantage of the structure of the Einstein equations, Klainerman 
and Rodnianski 9 were able to prove that, for CMC (constant mean curvature) foliations, smooth 
solutions exist beyond time ti if 

sup (||k||ioo(„^) + ||Vlogn||ioo(„^)) < 00. 
te[to,ti) ^ ^ 

This result was improved by Wang |12) who assumed an integral rather than a sup-norm condition: 

J^^ (l|k||L-(«,) + ||Vlogn||L=(H,)) dt < 00. 

Based on the results established in the present paper, especially our new inequality for the Bel- 
Robinson energy, we conjecture that the breakdown criterion for the Einstein equations can be 
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drastically weakened and a one-sided bound on nvr should be sufficient, that is, with the notation 
in Theorem \TM 

I sup(p(n7r)) dt < oo. 

J to 

In particular, this conjecture could be investigated in the case of CMC foliations. This objective 
seems to be realistic, based on the additional observation (made here by the authors) that the 
following a priori bounds hold (in terms of the Bel-Robinson energy, only): 

lk||L6(„^) + ||D7r||i2(„^) < 1. 
Wn-'h^^n,) + WV^Wmn,) + || V'(n-i)||i2(„^) < 1, 

IMm-Ht) + l|Vn||L2(„^) < 1. 

Finally, we refer to Grant and LeFloch [51 15] for the derivation of injcctivity radius estimates for 
Lorentzian manifolds that enjoy a one-sided bound, only. 
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